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Abstract 

We consider a class of simple quasi one-dimensional classically non-integrable 
systems which capture the essence of the periòdic orbit structure of general hy- 
perbolic nonintegrable dynamical systems. Their behavior is simple enough to 
allow a detailed investigation of both classical and quantum regimes. Despite 
their classical chaoticity, these systems exhibit a "nonintegrable analog" of 
the Einstein-Brillouin-Keller quantization formula which provides their spec- 
tra explicitly, state by state, by means of convergent periòdic orbit expansions. 
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I. INTRODUCTION 



Very few quantum systems can be solved explicitly. Among them are the Standard 
textbook examples, such as the harmònic oscillator or the hydrogen atom [|I| . In all of these 
cases the spectrum of the quantum system is obtained as an explicit analytical formula of 
the form u E n = where n is the quantum number of the system. This procedure fails 
already for some of the simplest quantum systems, which are still considered elementary 
textbook problems. An example is a quantum particle in a box with a step potential inside, 
as shown in Fig. 1. Even for the simple problem of Fig. 1, explicit analytical solutions of the 
form u E n = ..." are no longer available since the problem leads to a transcendental spectral 
equation. The recommended method of solution is either numerical or graphical fH-fl]. We 
recently found a way [|]-[(| of obtaining explicit analytical solutions of a wide class of problems 
such as the one shown in Fig. 1, thus obtaining an explicit analytical solution of textbook 
problems which until now were relegated to numerical or graphical solution techniques. Our 
methods are also a step forward in the mathematical theory of almost periòdic functions f?j 
since we obtain explicit formulae for the zeros of a wide class of almost periòdic functions. 
Furthermore, the classical dynamics of the quantum systems to be discussed in this paper is 
chaotic. Since it may well be true in general that the quantized versions of classically chaotic 
systems do not admit for the existence of quantum numbers (see, e.g., |8||5|] for a detailed 
discussion of this important point), our u E n = ..." spectral formulae, containing an explicit 
quantum number n, may come as a surprise. At this point we feel that it is important to 
stress that our results are not conjectures, approximations or merely formal identities. Our 
results are exact, explicit, convergent periòdic orbit expansions that can be cast into the 
form of mathematical theorems. We will publish the rigorous mathematical underpinnings 
of our results elsewhere [|ÏOfl . 



It is well known |]TTJ that the periòdic orbit theory leads to completely different ap- 
proaches for quantizing integrable and non-integrable dynamical systems. For integrable 



systems there is a simple procedure [11,12] which allows us to quantize the action vari- 
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ables individually for each degree of freedom. The situation is completely different for the 
chaotic case, where the periòdic orbit theory [D] allows us to evaluate only certain global 
characteristics of the spectrum, e.g. the density of states 

OO 1 oo 

p{E) =Y,&{E- Ei) « p{E) + - lm^T p (E) £ A v p {E) e^ E \ (1.1) 

j=l 71 p v=l 

typically only with semiclassical accuracy ||13|| . Here p(E) is the average density of states, 



S P (E), T p (E) and A P (E) are correspondingly the action, the period and the weight factor 
of the prime periòdic orbit labeled by p, and v is the repetition index. In this approach 
individual energy levels are obtained indirectly as the singularities of the sum ( |1 . 1|) . As for 
the idea of expressing them directly in terms of the periòdic orbits, M. V. Berry wrote in 



1991 |pH) : "...We do not know how, or even whether, the closed orbit sum generates the 
individual ós in the level density for chaotic systems. This is a serious - perhaps shocking - 
situation, because it means that we are ignorant of the mechanism of quantization". 

In the case of quantum graphs Berry's question can be answered defmitively. The periòdic 
orbit sums representing the spectral density of quantum graphs do provide the individual 



levels in the form of 5-spikes in flL·lp and only those [r5|-|IB|. In addition we showed recently 



that the answer to Berry's question can be taken one step forward: not only do periòdic 
orbit expansions for quantum graphs produce í-functions for the quantum states in the level 
density, but for certain classes of quantum graphs there also exist explicit convergent periòdic 
orbit expansions for individual energy levels. Because they provide explicit formulas for the 
energy levels of classically chaotic systems, these periòdic orbit expansions may be considered 
as "nonintegrable anàlogues" of the Einstein-Brillouin-Keller (EBK) quantization formula 



ïï| , p!2[| which applies to integrable systems. 

This paper is organized as follows. In Sect. II we briefly review the theory of quantum 
graphs and extend the theory by defming "dressed graphs", i.e. quantum graphs with 
arbitrary potentials on their bonds. In Sect. III we define an important class of dressed 
quantum graphs: regular quantum graphs. Based on a detailed study of their spectral 
properties in Sect. III we derive explicit analytical spectral formulae for regular quantum 



graphs in Sect. IV. In Sect. V we present a variety of regular quantum graphs illustrating 
use and convergence of the spectral formulae. In Sect. VI we summarize our results and 
conclude the paper. 

II. D YN AMIC AL NETWORKS 

Let us consider a particle moving on a quasi one-dimensional network of bonds and 
vertices. In the mathematical literature these networks are known as graphs. They were 
and still are the subject of intensive investigations in all areas of science ranging from 
mathematics over computer science to chemistry and physics. An example of a simple 
graph with five vertices and seven bonds is shown in Fig. 2. The particle scatters randomly 
at every vèrtex Vi along different bonds B^ which meet at that vèrtex. We shall assume 
that the graph contains a finite number of bonds and vertices (Nr and Ny correspondingly). 
The key assumption about the dynamics of the particle is that the turning points of any 
trajectory of a particle moving on the graph coincide with the vertices of the graph, and 
hence the shape of the trajectòries is determined uniquely by the geometry of the graph. 
The trajectòries of the particle are simply the joint sequences of graph bonds, which are 
easily described and enumerated. For instance, every trajectory can be represented by a 
sequence of Ng symbols, each one of which corresponds to a certain bond [HJ. Since the 
trajectòries correspond to various bond sequences, every trajectory is described by a code 
word consisting of Nb symbols. 

We "dress" the bonds Bij of the graph with potentials, Uij(x), which may affect the way 
a particle moves along the bonds. However, it is required that these dressings do not violate 
the geometry of the particle trajectòries, i.e. do not add turning points, other than the 
original vertices of the graph. This condition is required to hold at all energies. To comply 
with this requirement the bond potentials are allowed to depend on the energy E of the 
particle, i.e. £4, = Uij(x,E), such that E > Uij(x,E) is fulfilled for all E and all i, j. This 
will in fact lead to many additional simplifications which have a deep physical meaning in 
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the context of the semiclassical periòdic orbit theory [^-|6UÏ9|-pï| . 

The shapes of the trajectòries, and in particular of the periòdic orbits, becomes increas- 
ingly complicated as their lengths grow. This is what makes them similar to the genèric 
(dynamical) chaotic systems. In fact, the number of possible periòdic orbits increases expo- 
nentially with their lengths, (or, equivalently, the number of vèrtex scatterings) with a rate 
which depends only on the topology of the graph. Every graph T can be characterized by 
its topological entropy (global average rate of exponential proliferation of periòdic orbits) 

Ar = lim P.D 

where l characterizes the lengths of the periòdic orbits in terms of the lengths of their code 
words and #(/) is the total number of periòdic orbits of length < Z @. Since the phase space 



of the system is bounded, the dynamics of the particle is mixing fl6| . Hence the structure 
of the periòdic orbit set on dynamical networks closely imitates the behavior of the closed 
trajectòries of genèric chaotic systems |2^,[23|]. On the other hand, dynamical networks can 



be easily quantized |^-|6|,^7 HÏ9| , p4| , which makes them very convenient models for studying 
various aspects of quantum chaology. 

The details of the classical dynamics on graphs are discussed in numerous publications 
|Ï6| , |25[| . Below we investigate the quantum-mechanical description of these systems. In 
particular, we shall discuss their spectra in the context of the periòdic orbit theory. So let 
us now briefly outline some details of the graph quantization procedure, which will be used 
in the subsequent discussion. 

A quantum graph system is a quantum particle which moves on a one-dimensional net- 
work T dressed with the potentials Uij(x,E). Below we shall consider the case of scaling 



potentials discussed in P,pb[- 

Uij(E) = XijE, Xij = Xji, (2.2) 

where the X^s are constants. This choice of the dressing potentials allows us to avoid certain 
mathematical complications, which are irrelevant for the physical context of our discussion. 



For more details on scaling potentials and their relevance in the semiclassical periòdic orbit 
analysis see PJ^JT9|1. 

The Schròdinger equation for graphs with the potentials (p.2|) can be written as 



n&ijix) = plE^x), (2.3) 

where 

% = -ij- - Aij (2.4) 

is the generalized momentum operator and /3fj = 1 — Ay. The coordinate < x < Lij is 
measured along Bij from i to j, and = Lji is the length of the bond. The magnètic field 
vector potential = —Aji is assumed to be a constant, real matrix, and may be used as 
a tool for braking the time-reversal symmetry. 

Classically, the particle can travel along the bond if its energy is above the scaled 



potential height, E > Uij(E), (Ày < 1). In this case the solution of equation (|2.3| ) on the 



bond B^ is a combination of the free waves, 

gi( (3ij k-\-Aij )x 

ipij(x) = aij j= h kj == — , (2.5) 

where k = yE and the factors {f3ijk)~ 1 ^ 2 are introduced to separate the physically mean- 
ingful flux amplitudes from the coefficients and b^. In the opposite case of À v - > 1, the 
bond carries a linear combination of tunneling solutions. Due to the scaling assumption, 
there is no transition between these two CclSGS clS du function of E. From now on we shall 
assume that the energy E is kept above the maximal scaled potential height, 

\j<l, i,j = l,...,N v . (2.6) 

At every vèrtex Vi, the bond wave functions satisfy the boundary conditions, 

ipij(x = 0) = ifiCij, 

N v 

y2 ( ''j~'J C 'j(-''ï I x=o = -iK<Pi, (2.7) 
i=i 
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for all i,j = 1, ...,Nv- Here is the connectivity matrix of the graph, ipj is the value of 
the wave function at the vèrtex Vi, and the Àj's are free parameters of the problem, scaled 
as Àj = \®k (see Appendix). Note that the double-indexed scaling constants \j refer to 
the bonds, whereas the single-indexed constants A, refer to the scattering strengths at the 
vertices. We believe that this notation is natural and does not lead to confusion. 

The conditions fl2.7|) are consistent only for a discrete set of energy levels, E n = k\ which 
define the spectrum of the dressed quantum graph problem (|2.3| ) and ( [2. 71) . As shown in 
|6|JÏ6 Hr8l·p4| (see Appendix), using the scattering quantization approach f29|, one obtains 
the spectral equation for any quantum graph problem in the form 

A(Jfe) = det[l - S(k)\ = 0, (2.8) 

where S(k) is the finite unitary graph scattering matrix [TjJ. The indices that define the 
matrix elements Su of the matrix S correspond to the graph bonds. It is important that the 
bond Bj = Bij is considered to be different from the (geometrically identical) reversed bond 
Bp = Bji, so the bonds of the graph are directed P|-|5|,P^,P^[] . Hence the dimensionality of 
the scattering matrix is 2Nb X 2Ng. It is shown in the Appendix that S = TD(k), where T 
is a constant 2Nb X 2Nb unitary matrix and D is a diagonal unitary matrix, whose matrix 
elements are given by 

D IJ = 6 IJ e i ^ k+A ^ Ll , 7 = 1,..,2JV B . (2.9) 

Since A(k) is a complex function, it is convenient to define the spectrum via the zeros of its 
absolute value, 

|A(fc)| = e- i@o{k) A(k), (2.10) 

where Qo(k) is the complex phase of A(k). The logarithmic derivative of |A(fc)| produces a 
delta-peak for every one of its roots, 



1 d 00 

- Imlim — ln | det [1 — S (k + it)} \ = Y J 5{k- k n ), (2.11) 
7T ^o* - n=1 



which, by definition, is the density of the momentum states p{k) |J. On the other hand, 
using ( |2.10| ) and expanding the logarithm of the determinant (|2.8|) , the density of states can 
be written as 



Then it can be easily seen from the structure of the scattering matrix 5* that the 

matrix elements of its n-th power are defined on connected sequences of n bonds and the 
trace of S n generates terms defined on closed connected sequences of n bonds |||T7].^|.^] . 

These periòdic connected sequences of n bonds B^, can be viewed as the periòdic orb its 
traced by a classical point particle moving on the graph. Note that the phase of the exponent 
in ( p.9|) is exactly the action of a classical point particle's trajectory traversing the bond Bi, 

S x = í (p I k + A I )dx = (p I k + A I )L I . (2.13) 

Thus the semiclassical transition amplitudes e lSl between the vertices connected by the bond 
Bj determine the scattering matrix S(k). As a consequence fl-pi ^6l , ^9| > the "closed bond 



sequence expansion" (|2.12| ) can be written explicitly as a periòdic orbit expansion in terms 
of the phases ( |2.13|) , 



1 oo 

p(k) = p{k) + - Re £ S° p £ A; e^ k } (2. 14) 

71 V v=\ 



where S® is the fc-independent "action length" of the orbit p, 

S p = PijLijk = Sp k, (2.15) 
v 

and A p is its weight containing the constant factor A ^ L ^ . J3 ue to the scaling assumption 
(see Appendix) the weight factor A v is fc-independent. The first term in this expression 
corresponds to the average density of states of momentum p(k), 

m = If^fi, (2,6) 

while the periòdic orbit sum in Q2.14j) describes the fluctuations around the average. 
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The periòdic orbit expansion for the staircase function 

oo 

jV(fc) = £ 6(fc - fc„) (2.17) 

n=l 

can be obtained by direct integration of ( 2.1 1|) and (|2.14 ). We obtain 



N(k) = N(k) + N(k), (2.18) 

where the first term, 

N(k) = [ p(k') dk' + N(0) (2.19) 



o 



represents the average behavior of the staircase and 

i oo Au 

N(k) = Im - £ — e pk ( 2 - 2 °) 

n p u=l V 

describes zero-mean oscillations around the average. 

As discussed in the Introduction (see also |§-P,[Ï6],|25|), quantum graphs are chaotic in 
the classical limit. The classical scattering probabilities are obtained in the limit % — >• from 
the quantum mechanical transition amplitudes H-fl, (see Appendix). In the scaling case 
they are /c-independent and thus the quantum scattering amplitudes do not depend on h at 
all. They determine simultaneously the quantum and the classical scattering probabilities. 

III. REGULAR GRAPHS AND THEIR SPECTRA 

The spectral determinant is a polynomial of degree 2Ng of the matrix elements of S. It 
was shown in || that the total phase of this polynomial is 

e (Jfe) = ilmlndetS , (fc) = kS - 7r 7o , (3.1) 

where S = Lijfiij is the total action length of the graph T and 

1 Nv ÍX°\ 

7o = (N B + N v )/2 + -E arctan M- , (3.2) 

7T j=1 \Vi) 

where 



^2 CijPij. (3.3) 



Therefore the average density of states is a constant, 



?=^e„(*) = |. (3-4) 



and the average staircase ( |2.19| ) is 



N(k) = — k + N(Q). (3.5) 

7T 

The spectral equation |A(/c)| = can be written in the form 

N r 

cos (S k - 7770) = ^2 a i cos (Sik - Kji), (3.6) 

i=l 

where the freqüències Sj < S are combinations of the reduced classical actions S?- = ^Ly, 
and 70, 7$ are constants. The number N-p of terms in ( |3.6|) is bounded by N-p < 3^ ||. 

The frequency S'ò in Qo(k) is the largest frequency in the expansion (|3.6|) . While it is the 
only characteristic of the graph contained in the left-hand side of (|3.6|), the right-hand side 



^(k) =^2aiCos(Sik -TT^fi), (3.7) 
í=i 

contains the complete information about the graph system. We call $(&) the characteristic 
function of the graph. 

A graph V is called regular if its characteristic function $(/c) satisfies 

J2\ai\=a<l. (3.8) 
í=i 

For regular graphs the spectral equation ( |3.6|) can be solved formally 0-§J to yield the 

implicit equation of its eigenvalues, 

n ( arccos [$(&„)], for n + /i even, 

^ = [n + n + 7o] + { (3.9) 
D o >->o [ 7f — arccos[<&(/c n )], for n + fj, odd, 



where /i is a fixed integer, chosen such that k\ is the first positive solution of ( |3.6| ). The 
index n G N labels the roots of (|3.6|) in their natural sequence. 
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From the implícit form (|3.9|) it follows immediately that since the second term in ( |3.9|) 
is bounded by n/S , the deviations of solutions to this equation from the points 

L = + ^ + + (3.10) 

never exceeds 7t/Sq in absolute value for any n. Below it will turn out that the quantities 
k n are very important since they determine the root structure of (ft.6|) . 

The roots k n can be decomposed into an average part k n and a fluctuating part k n . From 
}.9|) we obtain 

k n = k n + k n , (3.11) 



where 



k n = ^[n + fi + l0 + l/2}, (3.12) 
>~>o 



K = ^ arccos [<£>(*;„)] - - } . (3.13) 



and 



So l /J 2. 

Note that the constant /i + 7o can be related to the initial value N(0) of the average staircase 
function ( |3.5| ). Consider the integral 

lim ^- / N(k')dk'. (3.14) 



n ^°° k n Jo 

The integration in Q3.14 ) can be easily performed due to the simple form ( 2.17 ) of the 



function N(k): 



1 rk n 1 J» 

^ N(k')dk f = n- ^J2 k ^ ( 3 - 15 ) 

K J ° K i= l 



since there are n roots to the left of k n . The fluctuations of both N(k) and k n around their 
averages have zero mean, so in the limit of n ^> 1 one can use N(k) and k n instead of N(k) 
and k n in ( p,15|) , and write 



1 rkn _ 1 JL _ 

^ N{k')dk' = n- ^J2 k ^ n>1 - ( 3 - 16 ) 

&n J ° K i=l 
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Using the explícit forms of k nj k n and N(k), we obtain 

JV( 0) + -(„ + /1 + 7o + l) = ___ ï (- + A + 7o + l). (3.17) 
Expanding the right-hand side and keeping terms up to order 1/n yields 

N(0) + ~(n + fi + 70 + 1) = n - (l - ^ + + Q + /i + 7 o + l) • (3.18) 
The terms proportional to n cancel out. Comparing the constants in (|3.18|) yields 



iV(0) = -(/i + 7 o + l)· (3.19) 
One can verify by direct substitution that 

N$ n ) = n, (3.20) 



which implies that the function (|3.10|) is the inverse of the average staircase (|3.5|) . One can 



also view the points k n as the intersection points of the staircase (|2.17| ) and its average Q3.5| ), 



N{k n ) = N{k n ) = n, (3.21) 
so at the points k n the fructuations N(k) of the spectral staircase vanish, 

1 OO AU 

N(k n ) = Im - E ~^e ivS > n = 0. (3.22) 

71 p v =l V 

Geometrically, ( |3.22|) means that the average staircase N(k) intersects every step of the 
staircase function N(k). Hence we call N(k) the piercing average. This is illustrated in 
Fig. 3 which shows the spectral staircase N(k) for the scaling step potential shown in Fig. 1 
and discussed in more detail in Sect. V, Example 1, below. We used the parameters À = 1/2 
and b = 0.3. Also shown is the average staircase N(k) for this case. It clearly pierces all the 
steps of N(k) providing an example of a system with a piercing average. 

Since $(/c) contains only freqüències smaller than S'ò, every open interval I n = (fc n -ij kn) 
contains only one root of ( |3.6| ), namely k n . Thus the k n play the role of separating points 
between adjacent roots. [P-|5|JT(]|] . Moreover, because of (|3.8|), the "allowed zones" R n C I n , 



where the roots k n can be found, narrow to 
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(Tï Tï \ 

— (n + fi + 70 + u) , — (n + n + j + 1 - u) , (3.23) 

where u = arccos(a)/S , . Correspondingly, there are forbidden regions F n , 

Fn = ( (n + H + 7o - u) , (n + n + 7o + 1 + uj) , (3.24) 

V>->0 ^0 / 

where roots of (OI) never appear. In the limit a — > 1 (u — > 0) the allowed zones i? n tend to 
occupy the whole root interval, R n — ► J n . 

IV. SPECTRAL FORMULAE 

Once the existence of separating points fc n has been established, it is possible to obtain 
an exact periòdic orbit expansion separately for every root of (|2.8|) . The derivation is based 
on the identity 



kp(k)dk. (4.1) 



n - 1 



Substituting the exact periòdic orbit expansion ( |2.14j ) for p(k) into (|4.1|), yields 

rk n Q n 1 rkn °° n <jr 

k n = I fc^dA; + - / fc Re £ 5 p ° A p e luS * k dk = —{n + /i + 7 o + 1/2) 

Jkn-l 7T 7T Jk n -i p u=1 Oq 

1 OO 41' 1 OO /If 

+ Re i £ J_ ^ fj ( e -^« - e lvS >^) . (4.2) 



Using (|3.22| ) we simplify ( |4.2|) to obtain 

O oo A v r , , n 

e íi/tjp(n+/i+70+l/2) g\ 



Tï 1 2 00 .A" 

fcn = ^(™ + ^ + 7o + 1/2) - -ImX) oo E -f sin 

where u; p = TïSp/S . The series expansion ( |4.3| ) for fc n is more than a formal identity. It is 
rigorously convergent, however it converges only condítionally, which means that the result of 
the summation depends essentially on how the summation was performed. Indeed, according 
to Riemann's well-known reordering theorem, one can obtain any result by rearranging the 
terms of a conditionally convergent series |ÏÏJ. Hence for proper convergence of ( |4.3| ) to the 
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exact roots of the spectral equation Q3.6D, we have to specify how the terms in Q4.3|) are to 
be summed. 



The mathematical details of the convergence properties of ( |4.3| ) are presented in |I0 
Here we mention the main result, which states that the terms in (|4.3|) have to be summed 
according to the length of the symbolic codes P,|H1 of the periòdic orbits, and not according 



to their action lengths. If (|4.3j ) is summed in this way it not only converges, but also 
converges to the exact roots k n of the spectral equation Q3.6|) . 

Hence the formula (|4.3| ) provides an explicit representation of the roots of the spectral 
equation ( |2.8|) in terms of the geomètric characteristics of the graph. According to ( [3.12|) , 



the first term in (4.3) is the average k n , and the following periòdic orbit sum is an explicit 
expression for the fluctuation of the root k n . This method is not limited to obtaining explicit 
analytical periòdic orbit expansions for k n . In fact, by using the identity 

f{kn) = L f{k)p{k)dk, (4.4) 

J k n -i 

we can obtain periòdic orbit expansions for any function of the eigenvalues / (k n ), for instance 
for the energy E = k 2 . 

In the simplest case where A° = 0, Ay = and lmA p = we have 

2 1 00 A v 1 
k n = ^n- - Y,^ T§ sin (o^p) sin^n). (4.5) 



"n ■ 



Note that in this case k^ n = —k Ti 

Both EBK theory as well as formula ( [4.3|) allow us to compute energy eigenvalues ex- 
plicitly. In this sense formula (|4.3j) may be regarded as an analog of the EBK quantization 



formula |TÏ1JÏ2|1 for a chaotic system. The complexity of this expansion, structurally similar 



to fll.ip , reflects the geometrical complexity of the periòdic orbit set for graph systems. 

Finally, for explicit calculations (see the following Sect. V) it remains to determine the 
explicit form of the expansion coefficients A p . For some simple graphs this was done in 
p~6| , p^9|1 . In the Appendix we solve the problem for general dressed graphs. We show that 
every passage of an orbit p from a bond Bij to B^i through a vèrtex Vi contributes a factor 
a ji,ij' ( a matrix element of the matrix T, see Appendix) to the weight A p of the orbit, 
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A P = II (4-6) 
where the product is taken over the sequence of the bonds traced by the orbit p. 

V. EXAMPLES 

In ( |3.8| ) we provided a definition of regular quantum graphs and discussed analytical 
properties of their spectra in (p.9|)-( |3T24l ). The discussion of regular quantum graphs culmi- 
nated in Sect. IV with the derivation of explicit spectral formulae for individual quantum 
states of regular quantum graphs. It is one thing to define regular quantum graphs, and it is 
quite another to show that regular quantum graphs actually exist. Examples 1-3, discussed 
below, provide specific examples of quantum graphs that are regular for all choices of their 
parameters. Examples 4 and 5 present quantum graphs that exhibit both regular and ir- 
regular regimes. Finally, examples 6 and 7 provide illustrations of a new class of quantum 
graphs, marginal quantum graphs, for which J2^=i \ a i\ — 1- Except for special choices of 
their dressing potentials these graphs can still be accommodated within the mathematical 
framework set up in Sections III and IV, and also admit an explicit representation of their 
spectra according to the spectral formulae derived in Sect. IV. 

Example 1. Scaling step potential inside of a box. Let us consider the case of a particle 

confmed to a box < x < 1 containing the scaling step potential (see Fig. 1 and Fig. 4a), 

f 0, for < x < b, 
U(x) = (5.1) 
[ A23-E, for b < x < 1. 

This is equivalent to a three- vèrtex linear chain graph (Fig. 4a'), with A2 = 0, = and 
Dirichlet boundary conditions at V\ and V3. This example is also discussed in [H-B, |I9 , 5T 
In this case the spectral equation Q3.6p can be written as ||,|5| 



sm 



HS°2i + S 2 3 )l = r sin \k(S° 21 - 5 2 ° 3 )1 , (5.2) 



where 

r = 

l + /?23 
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< 1 (5.3) 



is the reflection coeficient at the vèrtex V 2 . Hence the regularity condition ( |3.8|) is auto- 
matically satisfied and this graph is always regular. In Sect. IV we already discussed the 
convergence properties of (|4.3|), including the fact that a rigorous mathematical proof for the 



convergence of (|4.3|) exists |TÜ| . Here we present sòlid numerical evidence for the convergence 



of Q4.3D in the context of the scaling step-potential (|5.1|). As discussed in P~|5|,P^|, every 
periòdic orbit in the potential (|5.1|) can be described by a binary code word. Figure 5 shows 
the relative error e%} = \k$ — k n \/k n , n = 1,10,100, of the result k$ predicted by flOD 
compared with the numerically obtained exact result k n as a function of the binary code 
length l of the orbits used in the expansion (|4.3|) . We used 6 = 0.3, A = 1/2. Figure 5 also 
demonstrates that using all periòdic orbits up to binary code length l ~ 150, we obtain an 
accuracy on the order of 10~ 4 — 10~ 7 for the roots k n of ( |5.2|) . Although the convergence of 
the series is slow (according to Fig. 5 it is approximately on the order of ~ l/l 2 on average), 
one can obtain a fairly good estimate for the roots using all orbits of code length 20 and 
smaller. 

Example 2. Scaling 5-function in a box. This potential, shown in Fig. 4b, is again 
equivalent to a three-vertex quantum graph. This time, however, the potentials on the 
bonds are identically zero, whereas the vèrtex V 2 is dressed with a scaling í-function of 
strength A2 = X 2 k > 0. At the open ends we apply Dirichlet boundary conditions. In this 



case the spectral equation (|3.6| ) is 



cos 



where 



KS21 + S° 23 ) - VT70 = -\r\ cos k(S° 21 - S° 23 ) , (5.4) 



7o = 1 arcsin I = I (5.5) 



vr V v l • (A'M 

and the reflection coefficient r is given by 



A " (5.6) 



2i- \° 2 

Because of \r\ < 1, the characteristic function of ( |5.4| ) also satisfies the regularity condition 
( |3.8|) . Therefore the scaling 5 function in a box is another example of a regular quantum 
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graph. 

Example 3. Combined scaling step and scaling 5-potential in a box (Fig. 4c). This is 
equivalent to a three- vèrtex dressed linear graph (Fig. 4c'), with À2 = X 2 k > 0. The spectral 
equation (|3.6|) is 



cos 



k(S° 21 + S° 23 ) - 7770! = ai cos \k(S° 21 - S° 23 ) - 7T7il , (5.7) 



where 



7o 




71 = 1 arcsin (5.8) 



and the coefficient ai is 



Oi 



(fe ~ fe) 2 + (A^) 2 

So the characteristic function of (|5.7f) once again satisfies the regularity condition for any 
linear three- vèrtex graph with nontrivial bond potentials (/3 21 + /3| 3 7^ 0) ||. 

Quantum graphs which are regular for all of their parameter vàlues are quite exceptional. 
In general quantum graphs may have a regular regime for a certain range of the parameter 
vàlues or the regular regime may not exist at all. The following Example 4 illustrates this 
point. 

Example 4. Two scaling steps in a box (Fig. 4d). As an example of a graph which has 
both a regular and an irregular regime, let us consider a quantum particle in a box with 
two scaling steps (Fig. 4d), equivalent to the four- vèrtex linear graph shown in Fig. 4d'. 
Since there are no 5-functions present we have A2 = A3 = 0. We assume Dirichlet boundary 
conditions at the dead ends of this graph. In this case the spectral equation ( |3.6|) is given 
by 

sin(S , A;) = — r 2 sm(kSi) — r 2 r 3 sm(kS 2 ) + r 3 sin^S^), (5.10) 

where 
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C CO i QO i cO C CO i cO cO 

oo — o 2 \ -r o 23 i- o 34 , 01 — o 23 i- o 34 o 2 i, 

5*2 = S21 + 5*34 — S^g, ^3 = S 2 i + 5*23 ~ ^34 fó.ll) 

and r 2 , r 3 are the reflection coefficients 

Pl2 — @23 $23 — fíte / m m 

r 2 = o —o , ^3 = —5 (5.12) 

Pl2 + P23 P23 + P34 

at the corresponding vertices V*. For 

1^3 1 + |r 2 r 3 | + |r 2 | < 1, (5.13) 



the four- vèrtex linear graph (Fig. 4d', Fig. 6a) is regular. The regularity condition ( 5.13Q is 



fulfilled in a diamond-shaped region of (r 2 ,r 3 ) parameter space shown as the shaded area 
in Fig. 6b. The difference between the regular and the irregular regimes is clearly reflected 
in the staircase functions. Figure 7a shows the staircase function N(k) together with the 
average staircase N(k) in the regular regime for the parameter combination r\ = 0.2 and 
r 3 = 0.3. The piercing-average condition is clearly satisfied. Figure 7b shows the staircase 
function N(k) together with the average staircase N(k) in the irregular regime for the 
parameter combination r\ = 0.98 and r% = 0.99. In this case the piercing-average condition 
is clearly violated, consistent with the irregular nature of this regime. 

Example 5. Two scaling 5-functions in a box (Fig. 4e). This potential is equivalent 
with the four-vertex graph shown in Fig. 4e' with À 2 = X 2 k > 0, À3 = \®k > 0, and Dirichlet 
boundary conditions at the dead ends V\ and V4. In this case the spectral equation (|3.6|) is 
given by 

cos{kSo — 7770) = ai cos(kSi — irjx) + 02 cos^S^ — ^72) + a 3 cos(kSs — ^73), (5.14) 
where 

A2 A3 -^2^3 /r 1 r \ 

ai = a 2 = a 3 = = (5.15) 

V / 4 + (A^) 2 y/ï+W 2 v/[4 + (A^][4+(A^] 

and 
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1 . ( A°A°-4 \ 

70 = - arcsm , \ ; j , 

71 = iarcsin(a 2 ), 

72 = ^arcsin(ai), 

73 = i (5.16) 

The sum of the amplitudes in ( |5.15| ) ranges between and 3. Therefore this system has 
regular and irregular regimes. The regular regime corresponds to a finite area in the (Ai], A°) 
parameter space. All linear chain graphs with a finite number of vertices and Dirichlet 
boundary conditions at the two dead-end vertices at the beginning and at the end of the 
graph have a finite-measure regular regime and an irregular regime. This fact is proved in 

A new type of graphs are marginal quantum graphs. A marginal quantum graph is defined 

by 

£kl = i- ( 5 - 17 ) 

i=i 

For marginal quantum graphs, apart from a small set of "special" graphs, explicit spectral 
formulae still exist. Explicit examples are provi ded by circular graphs (see Example 6) and 
star graphs (see Example 7). 

Example 6. Scaling step potential in a box with periòdic boundary conditions. This 
system is identical with the two- vèrtex circular graph shown in Fig. 8. In the case of a 
circular graph a minor notational problem arises since starting from vèrtex V, say, vèrtex 
Vi may be reached directly via two different bonds. For the purposes of this example we 
solve the problem in the following way First we introduce a positive sense of rotation, 
i.e. mathematically positive or counterclockwise, for the circular graph in Fig. 8. Then 
we introduce the scaling coeflicient /3x2 referring the the bond that connects vèrtex V% with 
vèrtex V; traversing the graph in the mathematically positive sense. We introduce the scaling 
coeflicient /3 2 i which refers to the bond that connects vèrtex V2 with vèrtex V, again in the 
mathematically positive sense. We use the same notation for the two reduced actions S® 2 and 
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S 21 referring to the two different bonds (in the mathematically positive sense), respectively. 
This notation is not confusing here, since no magnètic field is switched on (Ay = 0). With 
this notation the spectral equation is given by 

cos(kS ) = ai + a 2 cos(kS 2 ), (5.18) 

where S Q = S° 2 + S^, & = S° 2 + S^, 

ai = — ; a , 2 (5.19) 



21, 



and 

fl2 = Ia7T& j ' (5 - 20) 

Note that a\ + a 2 = 1. Condition ( |5.17| ) is satisfied and the circular quantum graph with a 
scaling step potential is marginal. 

Although the strict inequality ( |3.8[ ) is violated, it is important to note that even in the 
marginal case the separating points k n , generically, are still not solutions to ( [5.181 ). This 
happens only for special parameter combinations, and therefore for special quantum graphs 
for which the equation 

(-l) n+ít+1 = ai + a 2 cos(k n S 2 ) (5.21) 

is exactly satisfied for some n. Since the sequence k n is countable, and (|5.21|) , in general, 
involves irrational requency ratios and irrational coefficients, equation ( |5.21| ) is only acci- 
dentially satisfied for some n for a measure zero set of graph parameters. Hence, in general, 
even for marginal quantum graphs, the points k n still serve as separating points and the 



roots of the spectral equation can still be obtained via the expansion ( [4.3[ ). 

Example 7. Star graph. Another example of a marginal quantum graph is provided by 
the star-graph shown in Fig. 9. We consider the case with three different scaling potentials 
on its three bonds and Dirichlet boundary conditions at the three dead ends. The spectral 
equation is given by 
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cos(So£;) = a\ cos(S\k) + a<i cos^/c) + 03 cos^S^k) 



(5.22) 



where 



C CO 1 CO 1 cO Q CO CO i cO 

J — ^14 t" °24 1" °34' °1 ~~ °14 °24 ' °34í 

C _ CO cO çO c _ oO 1 cO oO fe oq\ 

02 — o 14 o 24 o 34 , 03 — o 14 -t- o 24 o 34 ^Ü.ZOJ 



and 



^ _ /?14 — /?24 + @34 ^ _ — /?14 + 024 + /?34 ^ _ flli + /?24 — ^34 /g 
/^14 + 024 + /?34 /?14 + ^24 + /?34 /?14 + @24 + /?34 



It is straightforward to verify that no matter what the sign of the individual in ( |5.24|) 



Sf=i \ a %\ = 1- Therefore the condition (|5.17|) is satisfied and the star-graph shown in Fig. 9 
is another example of a marginal quantum graph. As discussed in the context of Example 
6, and with the exception of a set of measure zero of the star-graph's parameter space, 
the spectral expansion ( |4.3| ) is still vàlid and can be used to obtain every single one of the 
star-graph's eigenvalues individually, and independently of all the other eigenvalues. 

VI. SUMMARY, DISCUSSION AND CONCLUSION 

Exact periòdic orbit expansions for the global density of states are known for many 
chaotic systems |I5|,32|,33| . However, formula ( |4.3|) is the first example of an explicit expres- 



sion for the individual quantum mechanical levels, obtained as afunction of the level indexn, 
for a classically chaotic system. Additional explicit quantization formulas may be found for 
other quantum graph systems, or even for quantum systems unrelated to quantum graphs 
as long as two essential requirements are fulfilled. First an exact periòdic orbit expansion 
for the density of states must exist. Second it must be determined that one of the system's 
levels, fc*, is the only one in an interval < < k". Then one can always obtain the 
corresponding periòdic orbit expansion for 

rk'í 

K = I kp{k)dk, (6.1) 
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based on the periòdic orbit expansion for p(k). 

It is reasonable to expect that generically there exist separating points k[ and k" which 
separate every k* from its neighbors, so that fc* is the only root of the spectral equation in 
the interval Hence expansions like ( |4.3| ) generally do exist. However, knowing the 

positions of the separators k'^ and k" around a particular level fc*, does not help for finding 
the separators for all the other levels. So the most important task for obtaining a general 
expression for all the levels of a quantum chaotic system is finding a global function for the 
separating points similar to ( p.lOj ), which naturally enumerates the separators. Therefore, 
even though it might be possible to find, for some systems, the separators for a particular 
quantum level k*, and then to obtain a periòdic orbit expansion for it according to fl6.1|) , 
the expansion will work only for level k* and will not represent a formula which can be used 
for obtaining other levels. 

The problem of finding a global expression for the separating points as a function of 
their ordering index n is directly related to another well-known problem of spectral theory 
of differential operators, namely the problem of approximating the staircase function fl2.17| ) 
by a smooth average N(k). Indeed, suppose there exists a separating point k' n between every 
two roots of the spectral equation, i.e. a solution of the equation 

N(k') = N{k), (6.2) 

between every two roots of the spectral equation (similar to ( |3.20|) and (|3.21|) ). Since N(k) 
is a monotonic function, the separating points can then be found by inverting the equation 



k(N), (6.3) 



in which the value of the staircase function plays the role of the separator index k^. The 
equation ( |6.3j ) generalizes ( |3.10| ), which can be used in ( |6.1| ) to obtain the periòdic orbit 



expansions for all of the roots. 

The smooth curve defined by (|6.3| ) where N is considered to be a continuous variable, 
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intersects every stair of the spectral staircase ( |2.17|) . Unfortunately, finding a smooth func- 
tion which approximates the spectral staircase for a general differential operator with genèric 
boundary conditions is a rather complicated task. It was proven by Weyl in 1912 that one 
can approximate N(k) by the phase-space volume of the system in question, 

N(E) « JQ(E- H(x,p))¥^, (6.4) 

where D is the dimensionality of the phase space, however this average is certainly not 
guaranteed to satisfy the "piercing average" condition (|6.2j) . Since Weyl, this problem has 
been addressed by numerous researchers (see, e.g., f34j|), who succeeded in giving many 
improved estimates, for N(k), but none of them a priori satisfy ( p.2| ). 

The important feature of the regular quantum graph systems is that there exists a global 



piercing average (|3.10|) , which uniformly enumerates all the points separating one root from 
another, and hence it is possible to obtain formula (|4.3|) as a function of the index n. In 
other words, the index n in ( |4.3| ) is a quantum number, and so the expression ( |4.3j ) for the 
energy levels of a chaotic system in terms of classical periòdic orbits can be considered as a 
non- integrable analog of the EBK quantization scheme [ ]TT| , |T2"| . 

It should be mentioned that despite the existence of a quantum number n in ( |4.3| ), the 
actual dependence of the energy levels on the value of its quantum number is quite different 
from the simple EBK scheme for integrable systems. The expansion of the fluctuating part of 
the roots (|3.11|) involves an intricate, conditionally convergent series and is rather "chaotic" . 
The difference in complexity of the formulas ( |4.3|) and the EBK formula apparently reflects 
the complexity of the geometry of the periòdic orbits of the classically chaotic quantum 
graphs. 

Y.D. and R.B. gratefully acknowledge fmancial support by NSF grants PHY-9900730 
and PHY-9984075; Y.D. and R.V.J by NSF grant PHY-9900746. 
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VII. APPENDIX 



For the sake of completeness, we present below a simple derivation of the spectral deter- 
minant (j2.8|), starting from the boundary conditions at the vèrtex V{. 



ipij(%) \x=o= ViCij (7.1) 

and 



Í dx~ + Aij j ^ ij ( Xij } ' I=0= ^ iLf>i ' ( 7 ' 2 ) 



We would like to present the wave function, 

^(x) = — L= (<,,,< ;i J '' + l Hj e i3 « kx ) , (7.3) 
which satisfies these boundary conditions as a superposition of the partial waves 

Wr ( x j ) = s 3i' Tttt — + a w nrr ■ ( 7A "> 



scattering on the vertices of the graph. So 

N v iiPijk+Ai^Xj N v 

= E <VV#(^) = -p=e^ k -^ + E a i3' a ji,ij'i (7-5) 



with the appropriate weights a^i corresponding to the incoming flux on the bond Bjn towards 
the vèrtex Vi. Comparing this expression to fl7.3|) yields 

N v 

K = zJ a J·-·j'"'r· ( 7 - 6 ) 

i'=i 



Substituting (|7.5| ) into the boundary conditions (|7.1|) and (|7.2|) at the vèrtex Vi we get 
correspondingly 

li 7= ( ó Jf + = Vi c ij (7-7) 
j'=i yftijk 



and 



E Cijdij' y* [3ijk (8jj> - (Tji,ij>) = í\í<Pí, (7.8) 
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Substituting ( |7.7|) into ( |7.8| ) one gets 

Ny 



Ny 



üijl 



C {j J2 C ü a tj ,y%k Í5i f - crg,) = iXi 

l,f=i f=i \JPijk 

In case of the linear scaling, Àj = fcÀ°, this yields 

Ny Ny Ny f 

j'=i 1=1 j'=i y Hij 



Comparing the coeficients in front of a^v, we get 

r W „-\0 °3J' ÍA i 



1=1 



1=1 



Pij \ Pi 



or, after performing the summation over l, 



x „• \ o 

r r IR r v i ò\° — i 
^íj^ij'y Pij' ^ij 1 i j' lA i 



Pij \ Pi 



&ji,ij' 



where V\j, = £^ CuVPiíViiiï' ■ 

Multiplying both sides by CijPij and summing over j yields 



where Vi = J2j CijPij- Hence 



which can be used in ( |7.12[ ) to obtain 



r r ÍíT r Vl ~ lX ^ r HT a° 5 "' - iA » 9 

^ij^ij'V Pij' ^ij i ^Tïï^y'V Pi? i 



Vi + iX) 



&ji,ij' 



or 



2^J PijPiji 



We see that in the scaling case the matrix elements (Tjiiji of the vèrtex scatterinj; 
are /c-independent constants. 
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The matrix element crjijj has the meaning of the reflection coeficient from the vèrtex 
Vi along the bond and the elements <7ji t ij>, j 7^ j' are the transmission coefficients for 
transitions between different bonds. The equation ([7.6[) can be written as 



b = fa, (7.17) 

where 

In the symmetrical basis ipji (Ly — x) — tpij(x) one has 

e i(— Pijk+Aji)(Lij— x) e i(f3ijk+Aji){Lij-x) 

ipji (Lij - x) = ciji -== h bji j= = ipij(x), (7.19) 

so the coefficients and are related as 



CLji 



5 íj . e í (ftj íc + A ij) i ij ) 5 — a ^. e *(-ftj fc + A y)- L y_ (7.20) 



The coefficients and a^, (b^ and are considered different, so the bonds of the graph 
are "directed". 

Equations ( 7.20Q can be written in matrix form, 



a = PD(k)b, (7.21) 

where a and b are the 2Nb dimensional vectors of coefficients and D is a diagonal matrix in 
the 2Nb x 2Nb space of directed bonds, 

D ijm {k) = S ip S jq e^ k+A ^ , (7.22) 

and 

/ i Nb \ 

P=\ , (7.23) 
W B ) 

where 1n b is the iV^-dimensional unit matrix. The pairs of indices (ij), (pq), identifying 
the bonds of the graph T, play the role of the indices of the matrix D(k). 
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Equations (|7.21|) and (|7.17|) together result in 



a = S(k)a, 



(7.24) 



with the matrix S(k) (the total graph scattering matrix) given by 



S(k) = D(k)T, 



(7.25) 



where D = PDP and T = PT. 
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Figure Captions 

Fig. 1: Sketch of a step-potential in a box, a well-known textbook quantum problem. 
Fig. 2: Sample graph with five vertices and seven bonds. 

Fig. 3: The exact spectral staircase N(k) and its average N(k) for the scaling step-potential 
shown in Fig. 1 with b = 0.3 and À = 1/2. The average N(k) crosses every "stair" of N(k) 
(piercing average) at the equally spaced separating points k n . 

Fig. 4: Collection of potentials and their associated linear quantum graphs that serve as 
examples to illustrate the concept of regular quantum graphs. (a) Scaling step-potential 
in a box and its associated three- vèrtex linear graph (a'). (b) Scaling 5-function in a box 
and its corresponding three- vèrtex linear graph (b'). Combined scaling 5-function and step 
potential in a box (c) with its linear three- vèrtex quantum graph (c'). Two scaling steps (d) 
and two scaling 5-functions (e) in a box together with their associated four-vertex dressed 
linear quantum graphs (d') and (e'), respectively. 

Fig. 5: Comparison between the exact eigenvalues k n and the k n - vàlues computed via ( |4.3|) 
for the scaling step-potential shown in Fig. 1. Shown is the relative error e^' = \k^ — k n \/k n , 
n — 1, 10, 100, of the result k$ predicted by ( |4.3j ) compared to the numerically obtained 
exact result k n as a function of the binary code length l of the orbits used in the expansion 
O. We used b = 0.3, A = 1/2. 

Fig. 6: Four-vertex linear chain graph (a) and the corresponding space (r 2 , r 3 ) of reflection 
coefficients (b). The shaded region in the (r 2 ,r 3 ) space corresponds to the regular regime 
of the quantum graph shown in (a). This demonstrates that the subset of regular quantum 
graphs within the set of all four-vertex linear quantum graphs is non-empty and of finite 
measure. 

Fig. 7: The exact spectral staircase N(k) and its average N(k) for the regular r 2 = 0.2, 
r 3 = 0.3) (a) and the irregular r 2 = 0.98, r 3 = 0.99) (b) regimes of the four-vertex linear 
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graph shown in Fig. 6a. In the regular regime (a) the average staircase N(k) pierces every 
step of N(k). This is not the case in (b), characteristic of the irregular regime. 

Fig. 8: Two-vertex circular graph. In the mathematically positive sense, (3± 2 is the scaling 
coeficient of the bond connecting vèrtex V\ with vèrtex V 2 , foi is the scaling coefficient of 
the bond connecting V 2 with Vi. This labelling is possible only in the absence of a magnètic 
field (Aij = 0) where the sense of traversal of a bond is irrelevant. 

Fig. 9: Scaling star-graph with three bonds and four vertices. 

Fig. 10: Sketch of a piecewise constant potential ("Manhattan potential") (a) and its 
associated linear graph (b). 
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